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Abstract 

We review a lattice strategy how to non-perturbatively determine the coefficients in the HQET expansion of all com¬ 
ponents of the heavy-light axial and vector currents, including l/m^ -corrections. We also discuss recent preliminary 
results on the form factors parameterizing semi-leptonic B-decays at the leading order in 1 / m h. 
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1. Introduction 

B-meson decays present an excellent opportunity to 
test the Standard Model and to search for New Physics. 
Many decay channels are allowed. Their rates de¬ 
pend on CKM matrix elements and often have a sen¬ 
sitivity to heavy particles both within and outside the 
Standard Model. Such effects provide us with indirect 
tests of New Physics which are complementary to di¬ 
rect searches for new particles at very high energy. In 
this way, B-decays provide for stringent constraints on 
the structure and couplings of theories beyond the Stan¬ 
dard Model. Of course, to make this program work, the 
theory of these decays needs to be understood. Here the 
challenge is in the correct handling of QCD. 

The relevant intrinsic QCD scale is much smaller than 
the mass scale of the electroweak interactions or any 
New Physics. Therefore, the effect of beyond the Stan¬ 
dard Model interactions can accurately be described by 


a perturbatively treated effective Hamiltonian, consist¬ 
ing, e.g., of bilinear quark currents or local four-quark 
interactions. One then has to compute the hadronic ma¬ 
trix elements of this effective Hamiltonian with the basic 
theory given just by QCD (with N[ = 5 flavours). 

A first principles evaluation of the matrix elements is 
possible by lattice QCD, non-perturbatively in the QCD 
coupling but, as mentioned, perturbatively in the effec¬ 
tive Hamiltonian, which has contributions from the elec¬ 
troweak sector and, possibly, from New Physics. We do, 
however, still have to treat a hierarchy of scales 

L -1 « m n ^ 140MeV <$: m# « 5 GeV cT l . 

This forbids direct simulations on a lattice with spatial 
extent L and lattice spacing a on nowadays comput¬ 
ers. We therefore apply another effective theory, HQET, 
which provides us (at appropriate kinematics) with the 
expansion of QCD observables in inverse powers of m h. 
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The mass scale in the numerator of this expansion pa¬ 
rameter is the intrinsic QCD-scale of around 500 MeV 
as well as (in the rest frame of the B-meson) small spa¬ 
tial momenta. This theory, put on the lattice, only has to 
cover accurately the scales far below . 

HQET has low-energy constants (or HQET param¬ 
eters) which need to and can be determined non- 
perturbatively, as explained in m. Besides the HQET 
parameters in the action, there are a number of param¬ 
eters for each considered composite field which makes 
up the electro weak Hamiltonian. Extending the proce¬ 
dure discussed in (T| for the parameters in the action and 
the temporal component of the axial current, we review 
here a matching strategy for the full set of temporal and 
spatial components of the heavy-light vector and axial 
currents. 

These are of relevance for the determination of the 
CKM elements. For instance, V u b can be determined 
from the inclusive decays B —> X u £v (2l sect. 12] 

|y ub | = (4.41±0.15l°^)xl0- 3 , (1) 

or from various exclusive decays. From the exclusive 
semi-leptonic decay B —> nlv, one finds |2] sect. 12] 

| Vub | = (3.28 ± 0.29) x 1(T 3 , (2) 

while the purely leptonic channel, which is more diffi¬ 
cult in experiment, tends to give values of |V u b| that are 
significantly larger than in eq. For instance, using 
the branching ratio !B(B —> tv) = (1.14 ± 0.22) x 10 -4 
one obtains O sect. 12] 

|V ub | = (4.22 ± 0.42) X 10" 3 . (3) 

The discrepancies between these different determina¬ 
tions of | V u b | are not understood. We need both theoreti¬ 
cal and experimental improvements to find out their ori¬ 
gin. For exclusive decays, such as B —> tv or B —> nlv, 
it is known how to compute the required hadronic ma¬ 
trix elements by lattice computations. In this respect, 
the purely leptonic decay channel B —> tv, is the sim¬ 
plest case and the required decay constant /b is known 
precisely. We refer to Q and references therein for a 
summay. 

Semi-leptonic decays, such as B —> n£v or B^ —> KLv 
allow to extract the CKM parameter |V u b| from a com¬ 
bination of the experimental differential decay rate with 
a theoretical determination of the form factor f+(q 2 ). In 
particular lattice QCD can achieve the computation of 
f+(q 2 ). In principle, it is sufficient to do this at a single 
value of q 2 , the squared momentum transfer. However, 
in practice, experimental data is provided over a range 
(of bins) of q 2 , and one has to use some knowledge on 


the functional form of f+(q 2 ), e.g., in the form of the 
BCL paramterisation a Then, a theoretical prediction 
of f+(q 2 ) for a single q 2 is sufficient to extract |V u b|. In 
section 4, we investigate the feasibility of a precise lat¬ 
tice determination of the form factor in the continuum 
limit and at a fixed q 2 . Of course, eventually the preci¬ 
sion will be enhanced when the form factor is known at 
more than one q 2 . 


2. General matching strategy 


Following 0I61I3I8), we fix the 19 parameters in the 
HQET expansion of the action and vector/axial currents 
through matching equations, which we write in the form 

4>Q cd (L, m h , 0) = ®f QET (L, m h , a ). (4) 

In the above equations the O, (i = 1... 19) are finite- 
volume, renormalized, QCD quantities defined in the 
continuum 


® qcd (L, m h , 0) = lim ® qcd (L, m h , a) (5) 

1 0 1 


whereas the d> ; HQET are understood to be expanded up to 
a given order in 1 / m h (NEO in our setup) and computed 
in HQET at a finite lattice spacing. The 1 /m^-expansion 
of the observables O, results from expanding the action 
and composite fields 0 QCD (v), viz., 


0 HQET (m h ) = zJo stat + y\„0„ 


( 6 ) 


The coefficients in the expansion are chosen such that 
eq- © holds. In the r.h.s. of the above equation one has 
to include all linearly independent operators with mass 
dimension one higher than 0 QCD (or O stat ) transforming 
in the same way as 0 QCD under the common set of sym¬ 
metries of QCD and HQET. These are needed for the 
renormalization and 0(a) improvement of the effective 
theory, and in order to account in HQET for the mh- 
effects of QCD. Such effects are described by the mh- 
dependence of the parameters Z 0 and c n in eq. & Once 
the matching is performed at a given value of the heavy 
quark mass and fixed L = L\, the L-dependence in the 
difference ®^ CD (L, mh, 0) - ®^ QET (L, mh, a) quantifies 
higher-order (in 1/mh) corrections. These have been 
studied at tree-level in (8) for a particular class of match¬ 
ing conditions. Note that at tree-level and in the contin¬ 
uum limit, the L-dependence of dimensionless observ¬ 
ables <5>i appears only in the combination z = Lm h. 
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2.1. Definition of the matching observables 

Let us assume we want to determine the parameters 
appearing at NLO in the expansion of a current J. 

The matching observables can be constructed from 
suitable combinations of correlation functions, C/, 
which typically have a single insertion of J. Writing 
the HQET expansion of J as 


7 HQET = Z HQET I ystat + £ ^ l + 0 (l/ m g) , (7) 


then the correlation functions take the generic form 


£>HQET _ ^HQET^HQET^-mbare-^cj^stat _j_ 

^kinCy n + U; sp i n C^ Pm + ^ Cj n C l J ™J , (8) 


where all correlators on the r.h.s. are computed in the 
static approximation. The only appearance of the bare 
HQET quark mass, mbare, is in the factor e ~ mhaieXc , with 
xc given by the time distances of the heavy (static) 
quark fields entering in C^ QET . Apart from Z^ QET , all 
other (re)normalization factors contributing to C EQET 
are collected in Z^ QET . 

In the correlation functions on the r.h.s. of eq. ([8]), the 
leading-order term C^ tat has just one insertion of / stat (in¬ 
stead of /), while c^ n/spm differ from C^ tat by an extra 
insertion (summed over the entire space-time volume) 
of the 1/mh-terms 0\± n or 0 sp i n from the Lagrangian. 
The other next-to-leading contributions C l J m have an in¬ 
sertion of one of the higher-dimensional operators J n 
from the expansion in eq. ([7]). 

The observables O, are defined as combinations of 
such correlators, typically logarithms of ratios, and only 
one of the observables, say ®i, is left with an explicit 
dependence on mbare- 

By combining all HQET parameters into a vector 

CO — (mbare, ^kin, ^spin, Cj\ , • • • , hi Zj , . . .) , (9) 

the HQET expansion of the observables can be cast in 
the form 


0) f QET (L,M,a) = r h (L,a) + 

<p i i (L,a)a)j(M,a) + 0(l/ml), (10) 

where M is the Renormalization Group Invariant (RGI) 
heavy quark mass and the vector rj represents the con¬ 
tribution of the static terms C^ tat in the correlators in¬ 
volved. In general, the existence of a continuum limit 
can be guaranteed only for the combination of HQET 
quantities which appears on the r.h.s. of eq. ( fT0| ). Only 


at tree-level each individual term on the r.h.s. has a well 
defined continuum limit. 

In order to learn about the structure of the matrix cpi 
in eq. ( fT0| ), we group the HQET parameters into blocks 

(^bare), (^kin, ^spin), (O*, Zjf (C/', Z//), (C/", Zj"\ • • • , 

and assume that J is the current which is used in ®i. A 
suitable choice of the other observables then yields the 
following form of the matrix cpl: 

• All entries of the first column, except (p \, vanish. 

• In the first row, and are non-vanishing, due to 
the contributions from (o^ n and £U sp j n to <&\. In ad¬ 
dition, there may be non-zero <p J x with j from a sin¬ 
gle block, which corresponds to the current used in 
®i. In our case, this comes from the Ao,i-term (i.e., 
a term at NLO in the HQET expansion of the tem¬ 
poral component of the axial current, analogous to 
eq. (\3\ below), which enters the matching condi¬ 
tion for mbare- 


• The rest has a simple block structure, with non¬ 
zero blocks only in the second block-column (cor¬ 
responding to contributions from co^ n and uj sp i n ) 
and in the blocks on the diagonal (corresponding 
to the mixings in the last term of eq. ([8])). 

Schematically, we have the following block structure: 


<p\ 

* 

* 

0 

0 ) 

0 

* 

0 

0 

0 

0 

* 

* 

0 

0 

o 

* 

0 

* 

0 

0 

* 

0 

0 

* > 


That form can be preserved when an additional (effec¬ 
tive) operator J' is included in the matching. Therefore, 
the system can always be solved, and the HQET param¬ 
eters can be determined, by block-wise backward sub¬ 
stitution. 

Let us explicitly discuss the case of the (renormal¬ 
ized) heavy-light vector current in HQET. The temporal 
component is 


v 0 HQET M = T 


2 

C \x)+Y J Cy OJ Vo,i(x) 

i= 1 


, (ID 


with the static term 


y; tat W = (12) 

and two additional dimension-four contributions 

Vo,i(x) = (13) 

Vo,2<X> = ^hW, 


(14) 
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where all derivatives are defined from symmetric near¬ 
est neighbour differences and if/b and if/e denote static 
and light (relativistic) quark fields, respectively. In or¬ 
der to have access to a variety of different kinematical 
configurations, we will make extensive use here of gen¬ 
eralized periodic boundary conditions for the fermions 

i/f(x + Lk) = e Wk i//(x ), 1/f(x + Lk) = Hx)e~ Wk . (15) 

The phases Ok |9j|T0'|, sometimes referred to as twist¬ 
ing Ill] [12), can be employed to inject a momentum 
|0h - 0(:\/L in correlation functions mm. For exam¬ 
ple, it is clear that for 0 h = Of composite fields such 
as Vo, 2 M above can be associated with total derivative 
operators and, therefore, correlation functions at zero 
external momenta are not sensitive to their insertions. 

Considering now the spatial components of the vector 
current, the HQET expansions (see fl3l ) can be written 
as 


yf QET (x) = zJ QET 



with the following four extra terms: 

V*,i(*) = hO), 

Vk, 2 <x) = ^(x)l(Vf-^|)l/r h (x), 

Vk,3(x) = 4r { (x)^(V?+^f)yiy k i(rh(x), 

V kA (x) = ^(x)l(Vf + ^|)l/r h (x). 

The axial vector current components are listed in (8) . 
The classical values of the coefficients are 


C V 0 ,1 - Cv 0 ,2 - c v w 
and c Vk2 = c VtA ■■ 


C V M 

1 

m h 


1 

2 m h ’ 


(17) 


Concretely, as in EEIEI, it is advantageous to con¬ 
struct the matching observables using Schrodinger func¬ 
tional (SF) homogeneous boundary conditions at Vo = 0 
and Vo = t GHE3G6). Two categories of QCD corre¬ 
lation functions are formed from composite fields in the 
bulk, 0 < vq < T, and boundary quark fields £... £ 


• Boundary-to-boundary correlators 


~^JzY (^( u )ys£r(v) £V(y)75&(z)) . 

u,v,y,z 

K l (0(,Ob) = 

12 

~^eY Y {^( u )r^;(v)4(y)r^(z)), 

k u,v,y,z 

F Vo (Vo, Of, Of, Oh) = 

~^j6 Y (4( u )r5^(v)(y,)o(x)4(y)r5^(z)), 

u,v,y,z,x 

• Bulk-to-boundary correlators 
k Vo (*O,^,0b) = 

Y Y (( y i)oW4(y)r^(z)), 

k y,z 

ky(x 0 ,0 ( , 0 b ) = 

~^YY {( y i)'t( x )4(y)y^(z)}, 

k y,z 

(vo, Of, Ob) = 

-y Y {( y i)iW4(y)yi^(z)). 

yz 

where the label b refers to heavy relativistic quarks of 
mass close to the b-quark mass, and the subscript I 
stands for 0(a) improvement, as discussed in flOj 171 
for Wilson quarks. The correlation functions are illus¬ 
trated in m. where it is also explained how they relate 
to quantum mechanical Hilbert space matrix elements. 


E 

C 


0.55 

0.5 

0.45 

0.4 

0.35 

0.3 


( 0 ^ 02 ) = ( 0 . 5 , 1 . 0 ) 
(0 1 , 0 2 ) = (0.0, 1.0) 
( 0 ^ 02 ) = ( 0 . 0 , 0 . 5 ) 


I 


0 0.05 0.1 0.15 0.2 0.25 0.3 


1/z 


Figure 1: Tree-level continuum results for the parameter Wkin in 

j^HQET 


F i(0 ( ,0 b ) = 

Y {^( u )y5^b(v)4(y)y5^(z)>, 

u,v,y,z 

F{\0 e ,0e') = 


In order to determine the parameters a) i9 1 < i < 19, 
the above correlation functions as well as some re¬ 
lated to the axial current are combined to 19 observ¬ 
ables ii. First of all, ratios are formed such that 
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Figure 2: Tree-level continuum results for the parameter lnZ^ of the spatial components of the vector current, at static order versus 1/z (left) and 
at 0(1 1 nth) versus 1/z 2 (right). Q\ = 0 is clearly preferable, while already at 0 = (1,1,1) the corrections become very substantial. We also note 
that at tree-level and for 6 1 = 0 the dependence of the matching equations involving a single angle on any other parameter vanishes. Hence, the 
corresponding parameters are independent of whether or not using the ”EX” setup, and higher-order contributions are absent. 


the renormalization of the boundary quark fields drops 
out. Then, one searches for combinations and particular 
choices of 0/ such that the matrix ip has many zeros and 
is easily inverted^ Finally, one looks for choices where 
higher 0(1 /m^) terms are rather suppressed. Such a 
search has successfully been carried out at tree-level of 
perturbation theory (8), and we discuss some of the re¬ 
sults in the next section. 



Figure 3: Tree-level continuum results for the parameter cy kl (de¬ 
pending on two angles) of the spatial components of the vector cur¬ 
rent. Smaller values of ( 61 , 62 ) tend to yield smaller 1/z-corrections, 
and part of the higher-order 1/z-dependence (being smaller in the 
”EX” setup) is induced by other parameters. 


1 More generally, the matrix ip should have a good condition num¬ 

ber. 


3. Tree-level study of HQET parameters 

We now illustrate a few results of the tree-level com¬ 
putation (8J of the HQET parameters for the vector cur¬ 
rents that are needed for the semi-leptonic B-decays. 

After choosing a suitable set of observables O,, the 
HQET parameters are computed from eq. 0 inserting 
the HQET expansion eq. ©• At tree-level, each term 
in eq. © separately has a continuum limit. Since we 
are here not interested in the cutoff effects, the limit is 
taken term by term, and we then solve the matching 
equations 

Y J V J i -0Jj = ^ CD -rii 08) 

j 

at a = 0 for the unknown HQET parameters cOj. Due 
to the block structure of p and because of additional 
simplifications, which occur at tree-level, and with the 
choice of twist angles according to the strategy of (8), 
one can determine the parameters one after the other 
from the matching equations f]~8j ). 

A main motivation of the tree-level study in K9 was 
to investigate the effect of higher-order corrections (in 
1 /nth) in the system of matching equations and to test 
the choice of observables. 

At tree-level and in the static limit, the HQET pa¬ 
rameters at a = 0 must reproduce the exactly known 
classical values. After multiplying each parameter with 
the proper factors of m^ this limit is obtained for z = 
L —> 00 . The size of the higher-order contributions 
in 1 /mh can then be explored by looking at the rate with 
which the HQET parameters converge to their known 
classical values. One expects in particular 

= 1 + 0(1 /z), (19) 


c\ k l mh 
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^kin 

= 1+0(1 /z). 

(20) 

In Zyf ET 

= 0 + O(l/z), 

(21) 

In Zjf 

= 0 + O(l/z), 

(22) 

In z HQET 
mz % 

= 0 + 0(l/z 2 ). 

(23) 


Any 0(1 /z) deviations from the behaviour indicated in 
the equations above are to be understood as a sign of 
the 1 /mh-corrections neglected in the static approxima¬ 
tion of HQET. Once all 1 /rah-contributions are included 
in HQET and in its matching to QCD, as assumed here 
throughout, also the 0(1 /z) corrections to the HQET pa¬ 
rameters are fixed. Consequently, deviations from a lin¬ 
ear 1 /z-dependence signal higher-order corrections dis¬ 
regarded in the effective theory. They are also revealed 
by a non-vanishing dependence of the HQET parame¬ 
ters on the specific choices for the 0-angles on which 
the observables used in the matching step depend. 

Within our tree-level study, where the HQET parame¬ 
ters can be determined one after the other from eq. ( p~8| ), 
one can turn off the effects of 0(1 /m^) corrections in the 
previous steps by replacing all parameters determined in 
previous steps by their exact classical values. One then 
sees directly and only the effects of the new observ¬ 
ables considered in the present block. This procedure 
is dubbed ”EX”. We note that the difference between 
”EX” and the solution of the full system is dominated 
by the 1/mh-effects in (figure[l]). 

Figure [T] shows that has rather linear 1/z- 
corrections, which have little dependence on 0. The 
coefficient of l/z is around one, the standard magni¬ 
tude which is to be expected. However, for the investi¬ 
gated set of O f , the other HQET parameters show much 
smaller corrections, as we see below for a few examples 
and as is discussed in more detail in (8J. 

From figures [2}]4] one can conclude that 6\ = 0 or 
small is a good choice for observables depending on a 
single angle (such as lnZ^), whereas for the parame¬ 
ters via observables depending on two angles (such as 
cy k l ) the 0-dependence is rather moderate. The param¬ 
eters depending on three angles (such as cy k3 ) show 
a stronger influence by the 0-choices and sizable ef¬ 
fects of higher-order corrections in 1 /z, which often do 
not decrease substantially (or not at all) in the ”EX” 
case; thus, these effects appear to be a genuine prop¬ 
erty of them. Figure [4] also compares two alternative 
ways to determine lnZ^, either based on two- or three- 
point functions [8,18]. Owing to the visibly flatter In¬ 
dependence, the latter is obviously to be preferred, and 
In Zy^ ET and In Z*? QET behave similarly. For more de¬ 
tails we refer again to (8). The matching strategy of (8) 



Figure 4: Tree-level continuum results for the renormalization con¬ 
stant of the temporal component of the vector current, In Zy 0 , in the 
static approximation. We compare its definition via two-point correla¬ 
tion functions only (2pt) with our standard definition to fix this HQET 
parameter, which uses three-point functions (3pt). 


is also studied at one-loop order in perturbation theory, 
see CD for a preliminary report of results, which sup¬ 
port the findings of the tree-level study. 

Let us finally come back to the higher-order correc¬ 
tions of 6o>kin> which were found to be large compared to 
all others. They propagate to a large extent into higher- 
order corrections of the HQET parameters of A 0 and 14, 
which are the current components required for the cal¬ 
culation of the B( S )-meson decay constant and form fac¬ 
tors of semi-leptonic B-meson decays to light mesons. 
However, as also discussed in El, it is possible to reduce 
these higher-order effects in 1 /z by removing them com¬ 
pletely at tree-level via a straightforward modification 
of the matching observables. This will be implemented 
together with the non-perturbative matching computa¬ 
tion. Furthermore, there is the option to reduce the 
higher-order effects in by a generalized set of O,, 
with less zeros in p. While this option may reduce the 
impact of 1 /m ^-terms further, we emphasize once more 
that the proposal of ID does have 0(1) coefficients as 
expected generically. 

4. B -> K Form Factor 

In this section, we summarize the exploratory work 
of the ALPHA Collaboration to determine the form fac¬ 
tors for B s —> KTv decays non-perturbatively on the lat¬ 
tice with Nf = 2 sea quarks I2T1 l22l 20]. For this de¬ 
cay no experimental data is available yet. However, the 
presence of a heavier spectator s-quark renders the lat¬ 
tice computations technically simpler than for B —> niv. 
Thus, it offers a good basis to gain solid control on the 
systematic errors. 
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id 

T/a X (LIof 

a [fm] 

nijj; [MeV] 

m n L 

# meas. 

# target 

A5 

64 x 32 3 

0.0749(8) 

330 

4.0 

500 

500 

F6 

96 x 48 3 

0.0652(6) 

310 

5.0 

254 

500 

N6 

96 x 48 3 

0.0483(4) 

340 

4.0 

220 

500 


Table 1: List of ensembles used in (2D- 


The decay amplitude for B s —> K7V is proportional to 
IKbl times the hadronic matrix element parameterized 
as 


(K(pk)|V /j |B s ( jP b s )> = 


= Mv 2 ) 


Pb+Pk- 


2 2 

m B “ m K 


ct 


2 2 
0 wt-mt 

+ /o(g ) B % V , 


(24) 


where the two form factors fo(q 2 ) and f+(q 2 ) depend 
on q 2 = (p Bs ~ Pk) 2 , the invariant mass of the lepton 
pair. When neglecting the lepton masses, only f+(q 2 ) 
contributes to the decay rate. 

In order to keep the value of q 2 fixed in the con¬ 
tinuum extrapolation of the form factors, we employ 
flavour twisted boundary conditions, eq. for the 
s-quark. By choosing the vector of twist angles, 6 , of 
the s-quark, one can thus freely tune the momentum of 
the kaon. To remain in the rest frame of the B s -meson, 
the heavy quark is twisted by the same angle. 

We perform our measurements on gauge field en¬ 
sembles generated with Nf = 2 dynamical sea quarks 
within the CLS effort. They were simulated with non- 
perturbatively 0(a) improved Wilson fermions and the 
scale is set via /k [231, i.e., physical units are derived 
from /k = 155 MeV. All ensembles have m n L > 4 
where finite-volume effects can be neglected. In this 
work, we use the three CLS ensembles A5, F6 and N6, 
described in f23l . Our error estimates take into account 
correlations and autocorrelations m. 

In the exploratory study presented here, we restrict 
our analysis to the static approximation because the 
non-perturbative determination of the full set of HQET 
parameters is still in progress. Thus we set u^kin = 
£i> S pin = cy = 0. For the renormalization constants to 
be used at static order in eqs. ( fTT) and ( [16] ), we follow 
the lines of [13, 25 ] and write 

4 QET = C PS (MJA m )Z™ Gl (g 0 )Z$? A (g 0 ), 

Z» QET = Cv(M h /A m )Z™ al (g 0 ). 

The matching factors, Cps and Cy, are known at three 
loops in perturbation theory [26], and Z^rgi is known 


non-perturbatively GZ). The truncation of the theory at 
static order is expected to be a 10-20% effect. 

For an analysis in the static approximation, we only 
need to consider the two- and three-point correlation 
functions (while in the same measurement runs we 
also determine the various correlation functions that are 
needed for the extraction of the form factors in HQET 
at order 1 /m^) 

C K (v° - y°; p) = 

C*(*°-y 0 ; 0) = 
cl j(t K ,t B -,v) = 

The P qiq2 (x) are interpolating fields, such as 
0qi( x )T5*/%>(•*), f° r ^e mesons. Different levels of 
Gaussian smearing 1251 of the s-quark in the heavy- 
light meson, i.e., different trial wave functions, are 
indicated by the subscripts i or j. For any suitable 
smearing i, the ratio 


2 e - ip(x - y) (P su (x)P us (y)> , 
x,y 

(y)>, 


x,y 


Z '■*’ 


(x K — Xy) 


X (F SU (VK)V7(Xy)F5 S (VB)) 


*k; q 2 ) = 


^H, E K t K E B t B 


^C k (^k)C?(^b) 


e 2 e 2 


(25) 


will give the desired matrix element in the limit of large 
Euclidean times, = x^ - Xy and tn = Xy-x^: 


<K(Pk)|Vju|B s ( 0)> = lim (26) 


Alternatively, we can parameterize the fe- and In¬ 
dependence of the correlation functions as 

C K (f K ) = 2^ m) ) 2 e-^, 

m 

C E (r B ) = 

n 

m,n 


and determine {/c (m) , from a fit to C K , and 
{0f\ E^} from a combined fit to C 2 . and C®. 
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The matrix element of eq. (26) is then given by the fit 
parameter We include the first excited B s -state, 

but only the kaon ground state (i.e., we take m = 1 and 
n= 1,2). 


The ratio f^ t10 of eq. d25j) is shown as a function of t b 
at fixed t k = 20 in figure [o] For comparison, the value 
of (pp’ 1 '* resulting from the fit is indicated as a red band. 



Figure 5: The ratio // atl ° (blue points) and the fit result (pj}’ 1 - 1 (red 
band) for lattice N6, n = 0 and fixed 7 k = 20. 


In the rest frame of the B s -meson, the matrix elements 
have the form 


<K|V 0 |B S > = ^m 2 ), 

(K|V;|b s > = 


where the form factors are related to (/+, /o) by 

1 1 


/+ = 
/o = 


/?- + F — 

y2m Bs y2m B; 

yf2m ; 


Xm Bs -E K )f ± , (27) 


(m Bs - E K )fy 


K 

+ (E\ “ OT |)/-l- 


(28) 


The form factor /+ extracted from the fitted tp ^' 11 at 
different lattice spacings is shown in figure [6] When 
computing /+ in the static approximation of HQET, we 
are free to keep or drop terms of order 1 /m h in eq. 

To illustrate the corresponding 0(1/mh) ambiguity, we 
show in figure [6] (and [7]) two sets of data points: the 
upper one corresponds to using all terms in eq. 
the lower one to dropping the term proportional to f\\. 
This ambiguity will disappear, once all 0(1 /mh) terms 
of HQET are included. For both sets we show a constant 
continuum extrapolation and one linear in a 2 . The lat¬ 
ter has by far the larger error, and within this error it is 
consistent with the result of the constant extrapolation. 


.^.^ 



x Data 

1 ■ • Extrapolations 

-0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

2 /r 2 -10- 3 

a j fm 


Figure 6: Continuum extrapolation of our data for /+ at q 2 = 
21.23 GeV 2 . 


♦ HPQCD, ref. [13] 

■ • This work, preliminary 


♦ 

♦ 

-1-1-1-1-1-1- 

r - 

* . 

i-1-1- 


17 17.5 18 18.5 19 19.5 20 20.5 21 21.5 22 


q 2 / GeV 2 

Figure 7: Comparison of lattice QCD results at various values of q 2 . 

Figure [7] shows a comparison of our results from the 
linear continuum extrapolation of f+(q 2 ) with recent re¬ 
sults of HPQCD l29l (at their smallest a = 0.09 fm and 
m n = 320 MeV). 


5. Conclusions 

With the matching conditions between HQET at 
0(1/mh) and QCD reviewed here, it has been demon¬ 
strated that the solution of the extended matching prob¬ 
lem, including all components of the axial and vector 
currents, is feasible. Our lattice formulation basically 
amounts to solving a linear system of matching equa¬ 
tions with the 19 HQET parameters as unknowns. The 
matrix associated with this system has a simple block 
structure. The realization of this matching strategy em¬ 
ploys observables which are constructed from corre¬ 
lation functions in a finite volume with Schrodinger 
functional boundary conditions, and which would cor- 
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respond to mesonic correlation functions in large vol¬ 
ume. We exploit particular kinematic settings that are 
well suited for a straightforward implementation on the 
non-perturbative level by numerical simulations. 

Although our matching strategy (entirely expressed 
through renormalized quantities in QCD) is not neces¬ 
sarily restricted to Wilson fermions, we will use these 
in the future non-perturbative computation. According 
to our past experience with Wilson fermions 0[7], very 
fine lattice spacings and good numerical precision could 
be achieved. Moreover, with our tree-level study of the 
1/z-dependence of the HQET parameters we have il¬ 
lustrated that the linear behaviour already sets in for 
z > 8 in most cases, with slopes of 0(1). In fact, 
this finding is very reassuring, as it implies that for the 
non-perturbative framework followed by our collabora¬ 
tion with a matching volume of extent L « 0.5 fm 
and dimensionless heavy quark masses around z = 13 at 
the b-quark mass scale, higher-order corrections beyond 
the already included 0(1 /nth) ones are suppressed by a 
factor of about 10. 

The non-perturbative solution of the HQET matching 
problem is a crucial step towards a precision compu¬ 
tation of phenomenologically relevant hadronic matrix 
elements in B-physics. We presented here the current 
status of our effort to compute form factors for semi- 
leptonic B-decays using HQET on the lattice. In par¬ 
ticular, we considered the decay B s —> K7V and inves¬ 
tigated two different methods to extract the form factor 
f+(q 2 ) at a fixed value of q 2 = 21.23 GeV 2 from the 
plateau value of a suitable ratio of correlators or from 
a simultaneous fit to the functional form of the correla¬ 
tors. 

The authors of [20 ] performed a continuum extrapo¬ 
lation of the lattice data for f+(q 2 ) and found very small 
O (a 2 ) effects. While these preliminary results were still 
computed in the static approximation and an extrapola¬ 
tion to the physical pion mass has yet to be performed, 
the preliminary value of /+ is in agreement with the re¬ 
sults from other collaborations at the same quark mass. 

There is thus no sign of reduction of the difference of 
V ub extracted from semi-leptonic and inclusive/leptonic 
decays, but at the moment there are 10-20% uncertain¬ 
ties due to neglecting the 1 /rah-terms. These terms have 
to be included for a firm conclusion. Once the HQET 
parameters are known non-perturbatively, they will all 
be included in the analysis. Moreover, we plan to ex¬ 
tend the computation to B —> ntv decays, several values 
of q 2 , and Nf = 2 + 1 flavours of sea quarks. 
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